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The change-over from phase- to amplitude-fluctuation driven superconductivity is examined for 
a composite system of free electrons (Fermions with concentration tif) and localized electron-pairs 
(hard-core Bosons with concentration ng) as a function of doping - changing the total concentration 
of charge carriers (ntot = nF + 2ub)- The coupling together of these two subsystems via a charge 
exchange term induces electron pairing and ultimately superconductivity in the Fermionic subsys- 
tem. The difference in statistics of the two species of charge carriers has important consequences 
on the doping mechanism, showing an onset temperature T* of incoherent electron pairing in the 
Fermionic subsystem (manifest in form of a pseudogap), which steadily decreases with decreasing 
n to t- Below T* this electron pairing leads, in the normal phase, to electron-pair resonant states 
(Cooperons) with quasi-particle features which strongly depend on ntot- For high concentrations, 
where ns — 0.5, correlation effects between the hard-core Bosons lead to itinerant Cooperons hav- 
ing a heavy mass m p , but are long-lived. Upon reducing the concentration of charge carriers and 
consequently ns, the mass as well as the lifetime of those Cooperons is considerably reduced. As 
a result, for high values of «b, a superconducting state below T* sets in at a T c , being controlled 
by the phase stiffness = h 2 n v /m v of those Cooperons, where n p denotes their density. Upon 
reducing ntot, the phase stiffness steadily increases, and eventually exceeds the pairing energy ksT* . 
There, the Cooperons loose their well defined itinerant quasi-particle features and superconductivity 
gets controlled by amplitude fluctuations. The resulting phase diagram with doping is reminiscent 
of that of the phase fluctuation scenario for high T c superconductivity, except that in our scenario 
the determinant factors are the mass and the lifetime of the Cooperons rather than their density. 

PACS numbers: 74.20.-z, 74.25.-q, 74.20.Mn 



I INTRODUCTION 

The features which characterize classical low tempera- 
ture superconductors are the disappearance, above a cer- 
tain critical temperature T c , of a gap in the density of 
states (DOS) of the electrons, occurring simultaneously 
with the disappearance of the magnetic field expulsion 
(the Meissner effect) and standard Fermi liquid behavior 
in the normal phase above T c . A further characteristic 
is the practical impossibility to change significantly the 
value of T c upon changing the concentration of charge 
carriers, because of T c being largely determined by the 
density of states (DOS) at the Fermi level, which gener- 
ally is not expected to change much with doping. 

None of these features are observed in the high tem- 
perature superconductors (HTSC). The opening of a gap 
in the DOS occurs gradually, as a multitude of different 
experiments Q) show. This gap initially emerges in form 
of a pseudogap - a dip in the DOS - below a certain tem- 
perature T* , which, depending on doping, can be much 
above T c . Upon lowering the temperature and approach- 
ing T c , this pseudogap smoothly joins the superconduct- 
ing gap, as is evident from its angular variation near the 
Fermi surface. If the opening of the pseudogap and of the 
superconducting gap represent different physical manifes- 
tations of one and the same pairing mechanism, an inter- 
play between these manifestations of electron pairing in 
the two phases is to be expected and is in fact observed 



in form of remnant effects, such as: 

(i) A remnant of magnetic field expulsion is seen in 
form of a transient Meissner effect several tens of de- 
grees above T c , judging from the optical conductivity in 
the Tera Hertz regime^]. This points toward long-lived 
diamagnetic fluctuations, which have been attributed to 
the presence of long-lived diffusing vortices above T c [[|, 
as well as to phase uncorrelated diamagnetic regions 
which act as precursors to the true Meissner state below 
T C [Q. Experiments, invoking Andreev reflections to in- 
terpret the enhanced tunneling conductance in the pseu- 
dogap phase j^, |(|, point toward phase uncorrelated pair- 
ing above T c . 

(ii) Remnants of local electron-pairing in the c-axes 
optical response (orthogonal to the CuO?, planes) in 
the pseudogap phase, are seen in the superconducting 
phase. Similarly, the c-axes component of the electronic 
kinetic energy is getting reduced upon entering the super- 
conducting phase |1, provided the normal state exhibits 
pseudo gap features. The doping dependence of the c-axis 
penetration depth, being qualitatively similar to that of 
the basal plane, suggests that both are strongly influ- 
enced by the pseudogap features of the normal state 
Tunneling measurements have indicated remnants of the 
pseudogap which continue to coexist with the supercon- 
ducting gap in the superconducting phase. The latter 
disappears at T c , while the former remains . 

The question whether such findings favor or not a com- 
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mon origin of pseudogap and of the superconducting gap 
is presently still being debated ]i"l|, Q. 

The features in the HTSC, involving the interplay of 
the pseudogap phase and of the superconducting phase, 
are highly doping dependent. T* steadily decreases with 
increased hole doping, while T c shows an equally steady 
rise until the two approach each other. T c then bends 
over and follows the descent of T* upon further doping 
and approaching the optimal/overdoped regime. T c plot- 
ted as a function of the phase stiffness (determined by 
the square of the inverse penetration depth), involving 
the ratio of the density of superfluid carriers over their 
mass, shows a universal linear behavior fill . 

The main emphasis in this paper will be to analyse the 
doping dependence of T* and T c within a precursor pair- 
ing scenario. Very little about that doping dependence is 
known for such a scenario when based on single compo- 
nent systems, such as the negative U Hubbard model or 
the effective BCS Hamiltonians, extended to strong cou- 
pling. In such studies, doping is frequently introduced 
ad hoc, by assuming a doping dependent electron hop- 
ping or inter-electron attraction, in view of simulating a 
physics close to a Mott transition. In the present paper 
we shall examine such doping dependent effects without 
making any such ad hoc assumptions. As we shall see, 
considering a mixture of itinerant electrons (Fermions) 
and localized electron-pairs (hard-core Bosons) coupled 
together via a charge exchange term, is capable of re- 
producing the doping dependent features of T* and T c 
specified above. The essential new features introduced in 
this two-component scenario are the difference in statis- 
tics of the two components and the hard core features of 
the short range two-electron resonances (Cooperons) in 
the itinerant electron subsystem which result as a conse- 
quence of the charge exchange coupling. The variation 
of T c is dictated by correlation effects of the Cooperons. 

This two-component scenario has some similarity to 
the single-component scenarios with attractive inter- 
particle interaction mentioned above, in as far as it can 
be viewed as a two-fluid picture of existing preformed 
pairs and unpaired electrons Q. The electron-pair reso- 
nant states which result in such a scenario have certain 
features which are reminiscent of localized resonance im- 
purity states seen in the HTSC, and which arise when 
Cu atoms are substituted by non-magnetic atoms, like 
Yet, the electron-pair resonant states which we 
are considering here, have the essential potentiality of 
becoming itinerant and thus to lead to a superconduct- 
ing phase controlled by excitations of electron-pairs with 
finite momenta rather than pair breaking. 

In section II we discuss the interplay between phase 
and amplitude fluctuations in systems with precursor 
pairing, contained in the spectral properties of the 
Cooperon propagator. This permits us to make a con- 
nection with the phase fluctuation scenarios, which have 
been widely discussed in the literature. In section III we 



briefly outline the model and the Green's function formal- 
ism which we adopt in order to treat the hard-core nature 
of the resonant electron-pairs. In section IV we present 
the results for the temperature variation of the density 
of states over a wide doping regime and analyse the dop- 
ing variation of T* . In section V we explore the spectral 
properties of the two-particle excitations and compare 
the low and high density regime of the Cooperons as far 
as their quasi-particle properties are concerned. In sec- 
tion VI we discuss the thermodynamics of the pseudogap 
phenomenon in terms of the specific heat and entropy 
for different doping regimes. Finally in section VII we 
give a summary of our findings which are characteristic 
of precursor pairing systems, involving Fermionic as well 
as Bosonic charge carriers. 



II THE COOPERON WITHIN A PHASE 
FLUCTUATION SCENARIO 

Based on the different experimental results mentioned 
in the Introduction, which clearly indicate that the 
physics of those HTSC is not BCS like (over at least 
a wide regime of doping), it is widely agreed upon that 
the onset of the superconducting state in the underdoped 
regime in those HTSC should be controlled by phase 
rather than amplitude fluctuations Jl6|, |l7| , |l8| |. This 
supposes the existence of local superconducting droplets 
with a given phase, pre-existing above T c and evolv- 
ing into a macroscopic phase locked state upon enter- 
ing the superconducting state. Such a situation can be 
realized provided that the fluctuations of the phase of 
the macroscopic superconducting wave function are less 
costly in energy than the fluctuations of the amplitude 
of the electron-pairs, which describe pair breaking. To 
within a first approximation, this so-called phase fluctu- 
ation scenario is generally described within a hydrody- 
namic formulation of a spatially fluctuating phase 4> and 
its conjugate variable 5n s , which describes the spatial 
fluctuations of the superfluid particle density n s . The 
corresponding effective Hamiltonian is given by 

H= 1 -D^r + ^(Sn s r, (1) 

where D$ — fi 2 {n s /m s ) denotes the superfluid phase 
stiffness, x the compressibility, and m s their respective 
mass. The temperature which controls the phase order 
of such a system is given by ksT^ ~ D^a, with a being 
either given by the coherence length £ for 3D systems or 
by the interlayer distance in layered compounds such as 
HTSC. 

This phase fluctuation scenario is frequently discussed 
in conjunction with the so called BCS-BEC cross-over 
phenomenon p9| , where, as a function of the strength 
of the inter-particle attraction, one passes from a BCS 
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state at weak attraction to a Bose Einstein condensa- 
tion (BEC) of tightly bound electron-pairs in the limit of 
strong attraction. The physics for that has been widely 
studied on the basis of effective BCS and negative U Hub- 
bard Hamiltonians, aiming to treat the single-particle 
and two-particle features on the same footing. The pseu- 
dogap in such a scenario arises from short range electron- 
pairing, correlated over a finite time scale, comparable to 
the energy scale of the zero temperature superconduct- 
ing gap. According to a general theorem (due to Bo- 
goliubov) and based on the singular behavior of the oc- 
cupation number of electron-pairs with small momenta 
which signal bound states, such electron pairing ought to 
survive below T e |po|| . A possible experimental verifica- 
tion for that might come from the socalled peak- dip-hump 
feature in ARPES|pl|, which shows a spectral behavior 
upon entering the superconducting phase, where a sharp 
peak (related to superconducting correlations of the low 
energy excitations) emerges out of the broad incoherent 
background, characterized by a broad hump and repre- 
senting remnants of the pseudogap phase. 

An attempt to formulate the problem of local pair- 
ing as a prerequisite of superconductivity was made a 
long time ago by generalizing the mean field BCS formal- 
ism such as to cover the regime above as well as below 
T c [^2). Instead of using the order parameter, the prop- 
agator for the electron-pairs, also called Cooperons, is 
introduced and treated on the same level as the single 
electron propagator. Within such a formalism, and on 
a quite general basis, the pseudogap phase of the HTSC 
has been examined p3[, invoking the mutual feedback ef- 
fect between the single- and the two-particle properties 
via the introduction of some effective gap above T c . 

A different procedure was followed by proposing a spe- 
cific structure of the Cooperon propagator |24| of the form 

C(r, t) = C a exp(-r/r ) + C^e^^e'^}, (2) 

separating amplitude from phase correlations in an ad- 
ditive way. The first term, representing a rather rapidly 
decreasing function with r, describes local pair ampli- 
tudes which are treated in a time independent fashion. 
The second term describes the phase correlations, which, 
in the low frequency limit, are approximately expressed 
by C^exp^r/^T)) where £(T) denotes the tempera- 
ture dependent coherence length. Attributing that latter 
contribution of the Cooperon propagator to a 2D XY 
physics above the Kosterlitz-Thouless critical tempera- 
ture Tkt: establishes a link with the phase fluctuation 
scenario. The doping behavior of such precursor systems 
as the HTSC is then monitored by parameterizing the rel- 
ative weight of the coefficients C Q ,C^ together with the 
relative spatial extent of the two contributions of C(r, t) 
such that it describes: 

(i) an underdoped regime, characterized by primarily 
phase fluctuation controlled onset of superconductivity 
with a large temperature regime for the pseudogap phase, 



(ii) an optimal/overdoped regime, characterized by a 
gradual disappearance, upon increased doping, of the 
pseudogap phase and a superconducting phase, con- 
trolled by amplitude correlations. 

Considering the origin of the pseudogap as being due 
to pair fluctuations, also called precursor pairing (not 
to be confused with preformed pairs), the characteristic 
temperature T* below which the opening of the pseu- 
dogap occurs for such single-component scenarios with 
inter-particle attraction, scales with the strength of that 
interaction^, f28|. However, as far as the con- 

centration dependence of T c and T* is concerned within 
such scenarios, it invariably shows that both T* and T c 
follow the same trend p9|, |30| upon varying the number 
of charge carriers. This is clearly the opposite to what 
is found in the HTSC. Potential dimensionality changes, 
linked to the change-over from underdoped to overdoped 
materials, can not remedy this situation either, since as 
far as T* is concerned, it is determined by essentially 
local physics and thus independent on any dimensional 
aspects. As far as T c is concerned, even upon assuming a 
cross-over, with reduced doping, to a Kosterlitz-Thouless 
critical temperature behavior in the underdoped regime, 
Tkt would still follow the same doping trend as T*. It 
thus seems likely that correlation effects are indispens- 
able in determining the doping dependence of T c versus 
T* in such precursor scenarios. 



Ill THE MODEL AND THE TECHNIQUE 
EMPLOYED 

The scenario of a mixture of itinerant Fermions 
(band electrons) and localized hard-core Bosons (bound 
electron-pairs) will be described on the basis of the so- 
called Boson- Fermion model (BFM). This model presents 
a paradigm for interacting electron systems where two- 
particle resonant states are expected to occur due to the 
interaction of the electrons with certain Bosonic modes 
and where such two-particle resonant states act as pre- 
cursor to a transition into a superconducting state. The 
underlying physics behind this model fyj, as it was ini- 
tially conceived H^], is that of electrons strongly coupled 
to local phonons, which act as such Bosonic modes. This 
results in self-trapped entities, comprising the charge car- 
riers and the surrounding clouds of Bosonic excitations, 
in form of resonant pair states inside a system of itiner- 
ant electrons. Such a BFM scenario is not in any way re- 
stricted to electronic systems undergoing a superconduct- 
ing transition but ought equally well apply to electron- 
hole pairing in semi-conductors |33| , and low density nu- 
clear matter with isospin singlet pairing |34|. Moreover, it 
was in a similar spirit that such a Boson- Fermion mixture 
scenario has been derived recently for (i) the Hubbard 
model with intermediate repulsive coupling j35|, (ii) the 
exchange interaction between spinon singlets of resonat- 



ing valence bond (RVB) electron-pairs and holons|3(| an d 
(ii) for entangled atoms in squeezed states in molecular 
Bose Einstein condensates in traps j37|. More generally, 
the BFM has been employed in attempts to bosonize an 
intrinsically Fermionic systemj38|. 

The model Hamiltonian describing the Boson Fermion 
scenario is given by 

Ho = (D - p) cf^ + (A s - 2p) $>? + \) 



+ 



(3) 



where the localized hard-core Bosons are represented by 
pseudo spin-1/2 operators [pf , p~ , pf] and the itinerant 
electrons by [ci a ,cf a ]. v denotes the strength of the on- 
site hybridization between the two types of charge car- 
riers and t the hopping integral for the itinerant elec- 
trons with a band half width D. The full band width 2D 
will be used as the energy unit through out this paper. 
The energy level of the localized Bosons is given by As 
and the chemical potential p is chosen to be common to 
both species, such as to ensure total charge conservation, 
n-tot — Tip + 2ub- np and ub denote the occupation 
number per site of the electrons (including up and down 
spin states) and of the hard core-Bosons. 

The opening of the pseudogap and its temperature de- 
pendence on the basis of this model, for a fixed con- 
centration and upon neglecting the hard-core nature of 
the Bosons has been studied previously J39|, [II], pi] . In a 
study, based on the dynamical mean field approach Q, 
the hard-core nature of the Bosons could be taken into 
account. But then, the itinerancy of the Bosons could 
not be treated within such a scheme which restricted this 
study to purely amplitude fluctuations driven supercon- 
ductivity. In the present work we shall account for both, 
the hard-core nature of the Bosons as well as their po- 
tentiality becoming itinerant and shall study the pseudo- 
gap characteristics as a function of total carrier concen- 
tration. The present study follows closely the previous 
self-consistent diagrammatic approach [|9), but general- 
izes it such as to take into account the hard-core na- 
ture of the Bosons. We adopt for that purpose the dia- 
grammatic technique which had been developed for spin 
systems 1 43 , 44f (with their convention [p + , p~] = p z ) and 
for which it was shown that the usual Wick theorem had 
to be generalized to 
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FIG. 1: Diagrammatic representation of the Fermion (a), Bo- 
son (b) and the vertex correlation functions (c+d). 
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(5) 



Recursively applying this procedure of the modified 
Wick theorem, the remaining multi-spin correlation func- 
tions are transformed step by step into a sum of prod- 
ucts of K®j multiplied with thermal averages of the type 
(pf....p%)o (evaluated with respect to the unperturbed 
Hamiltonian: (As — 2/i) Y2i(Pi + f))- These thermal av- 
erages are expressed in terms of a set of cumulants, the 
first few of which are given by: 



(Pt(r)) 

(pf(r)pl(r)pl(r)> 
+b"6 1 , 2 6 2 ,3 



1 

b'S U2 



77. B 



b' 



= b 3 + bb'(S h2 + 62,3+63,1) 

, b 3 + 366' + 6" = -6 
4 



(6) 



Keeping only the first two cumulants 6 and 6' gives rise to 
the set of diagrams, illustrated in Fig. 1, and describe the 
set of self-consistent equations determining the Fermion 
(Fig. la) and Boson Green's functions (Fig. lb) G(k,oj n ) 
and K(q,u> m ). The vertex depicted by the full square 
in those figures is made up of two contributions: one 
arising from the cumulant 6 and depicted by a full circle 
and one arising from the cumulant 6', illustrated by the 
second contribution to this vertex. Both those vertex 
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contributions have to be determined selfconsistently. For 
the contribution arising from the cumulant b it simply 
is b = —1/2 + For the contribution arising from 
the second cumulant, given by b'V(ui) the selfconsistent 
equation corresponding to the set of diagrams illustrated 
in Fig. Id has to be solved. To within this approximation 
of cumulants, this leads to the following set of equations: 



G(k,u n ) = 
K(q,u> m ) = 



KiAr,r') = (T[p-(r)p+(r')]) 



N/3 



q,m 



iu) n - e k - E(k,w n ) ' 

b + VV{u m ) 

iujm - E - (b + 6'V(w m ))II(g, uo m ) ' 



G°(k,u n ) = — ; K°(u m ) 



V(uJ m ) = V (ijJ m ) + K (u m )^^Il 2 (q,uj m )K(q,uj m ), 



1 \ - U.(q,u} m ) 



N ^ ioj m - Eq 



(7) 



with Eq = Ab — 2fi and the bare electron dispersion 
£k = D(l - E(r^r 3 ) e ik < r '- r J>) - p. The self energies 
for the Fermions and hard-core Bosons are given by: 



where uj n = n(2n + l)/j3 and oj m = ir2m/(3 denote the 
Matsubara frequencies for Fermions and Bosons respec- 
tively, n and m running over all integers from — oo to +oo. 
The expressions for the self energies for the Fermions, 
S(k, Lo n ) and for the hard-core Bosons, II(q, u> m ) differ 
from the standard ones for ordinary Bosons by a change 
in sign due to the Wick theorem for hard-core Bosons. 
The effect of the cumulants however corrects this sign 
change in the end because of the presence of the factor 
b in the numerator of the Bose Green's function. Fixing 
n to t, Eqs. (7) are solved numerically on the Matsub- 
ara axes and the resulting Green's functions and selfen- 
ergy functions are then analytically continued onto the 
real frequency axes, via the usual Pade type procedure. 
The Green's functions G and K are linked to the oc- 
cupation numbers n£ = (2/(3) E w e~ tUnQ G(k, u> n ) and 



S(k, u) n ) = --gjy X! _ k ' Wm _ CJ «)^( ( 1' w ™) Bosons respectively, with n F ,B = j? J2 



A'(q, Lu m ) for the Fermions and 



F,B 



2 

II(q,w m ) = — ^2 G(q - k, oj m - u> n )G(k, u n ) (8) 



This set of equations represents a generalization of 
the usual self-consistent RPA equations for this BFM- 
problemj3^, ^| when restricting oneself to the lowest or- 
der approximation involving only the cumulant b. The 
contributions arising from higher order comulants bring 
in frequency dependent vertex corrections V(u>). Qual- 
itatively, the hard-core nature of the Bosons is already 
contained in the lowest approximation, due to the ap- 
pearance of the factor b in the expression for the hard- 
core Bose Green's function. Given the considerable com- 
plexity in solving these equations we have, for a restricted 
set of values for the temperature and the total particle 
concentration, compared the results which arise from the 
full set of equations (0, ||) with those arising form the low- 
est order approximation, i.e., keeping only the cumulant 
b. The results being qualitatively the same, and given 
our aim to present only very robust qualitative features 
of the physics we want to discuss here, we report in the 
following exclusively the results based on this lowest or- 
der approximation. 

The Green's functions for the Fermions and the hard- 
core Bosons are defined by 



Gi tj (T,T') 



-(TMT)ct(r')]} 

e »k(ri-rj)-tu, 

N/3 ^ 

k,n 



>(t-t') 



G(k,w„), 



k "k 



As a typical example for the present study we choose 
the energy level of the hard-core Bosons to lie in the 
center of the band of itinerant electrons (Ab — 1) and 
assume a small value of the exchange coupling constant 
(v 2 = 0.02). Requiring the chemical potential to lie 
slightly below the Bosonic level (// < 0.5), assures us that 
upon changing the total number of charge carriers ntot 
from 2 to 1, we recover a situation which, as far as the 
density of electrons is concerned, mimics the situation en- 
countered in HTSC over a wide doping regime. We shall 
for that reason adopt the terminology, widely used in 
connection with studies on the HTSC's, and refer to the 
doping regime 2 > n to t > no as the underdoped regime 
and ntot < no as the optimal/overdoped regime, with 
n,Q ~ 1.1 for our choice of parameters. We furthermore 
restrict the momentum summations over a ID Brillouin 
zone with 200 fc-points. This is justified, since we are 
interested here in only very general features of the pseu- 
dogap phenomenon. If our results, have any bearing on 
the physics of HTSC, they should apply to regions in mo- 
mentum space where the pseudogap phenomenon is most 
pronounced i.e., near the so-called hot spots around the 
M points in the Brillouin zone of the basal plane. There 
they could describe d wave symmetry pseudogap behav- 
ior along one of it's lobs along a direction [0, 0] — [0, ir] and 
their equivalents, traversing the M points. Our results 
could thus possibly be compared with ARPES spectra 
for wave vectors along such directions in k space, as well 
to transport measurements along the same directions. 



6 




FIG. 2: Comparison of the variation of the number of normal 
to hard-core Bosons as a function of the chemical potential 
for several temperatures. 



IV THE DOPING AND TEMPERATURE 
DEPENDENCE OF THE PSEUDO-GAP 

We present in the following the results of the solution 
of the set of self-consistent equations (0, ||), together with 
a self-consistent determination of the lowest cumulant 6, 
Eq. (6). This permits us to determine the electronic 
DOS as a function of temperature for different "doping 
rates", given by n to t, or alternatively, by the depletion of 
the Fermi sea away from half filling, given by (1 — rip). 
In this BFM scenario, doping influences both the number 
of itinerant electrons np as well as the number of bound 
electron-pairs uq. This is not an unrealistic premise as 
far as the HTSC are concerned, since it has been exper- 
imentally established that doping does not occur exclu- 
sively in the CuOi planes but involves also the dielectric 
layers between them. This is born out by XPS studies 
which permit to determine the relative change with dop- 
ing of the population of Cu + versus Cu ++ ions [H. Fur- 
ther indications that doping occurs in a multi-component 
system comes from measurements of the size of the Fermi 
surface volume ^6| which show that the universal curve 
for T c as a function of doping is shifted downwards in dop- 
ing as compared to its dependence on the chemical dop- 
ing rate. And finally, site dependent XAFS studies Ejj] 
show that in order for the superconducting phase to ma- 
terialize, doping must necessarily involve holes located 
outside the metallic CuC>2 planes. Further evidence for 



the existence of two species of different charge carriers, 
itinerant ones (giving rise to a Drude peak) and localized 
ones (giving rise to a peak in the far infrared regime), 
comes from reflectivity measurements [fl8| . 

Given our choice of the Boson level falling in the mid- 
dle of the band of itinerant electrons, fixes the Fermi level 
such that we have the situation of a half filled band for 
n to t = 2. Upon hole doping we move the chemical po- 
tential downwards from its value at D, which introduces 
holes in the electronic subsystem and at the same time 
diminishes ns- This variation of as function of p 
is illustrated in Fig. || for different temperatures. The 
bound electron-pairs, being hard-core Bosons lead to a 
fully symmetric situation for particle and hole doping for 
this choice of parameters. We also illustrate in this fig- 
ure the variation of tib with p for ordinary Bosons, which 
significantly differs from that of hard-core Bosons. 

In order to study the evolution of the pseudo-gap as 
a function of temperature and doping, we evaluate the 
spectral function of the single-particle Fermionic Green's 
function Ap(k, ui) = 2ImG(k,ioj n — » uj + id)) which, af- 
ter integrating over all wave- vectors in the Brillouin zone 
gives us the DOS, p{u>). In Figs. || we plot the evolu- 
tion of the pseudo-gap near the chemical potential (cor- 
responding to uj = 0) as a function of temperature and 
for several doping concentrations n to t =— 1.66, 1.20 and 
0.97. In order to determine T* as a function of tib and 
ntot respectively, we illustrate in Fig. ^ the minimum of 
the dip in the DOS given by p m in{T), as a function of 
tib for different temperatures. We then consider the rel- 
ative values of this depletion of the DOS, determined by 
pminiT) I p m in(oo) and cut these functions by horizontal 
lines, lying 4% below the saturation values of p m in{oo). 
The crossing points determine the values of T* for any 
particular n to t and its corresponding value of ns- J 1 *, 
representing a cross-over rather than a phase transition, 
thus corresponds to that temperature where the devia- 
tion from the high temperature saturated DOS close to 
the Fermi energy is reduced by an arbitrary but small 
amount, chosen here as 4%. In Fig. [| we illustrate the 
variation of T* as a function of ntot as well as of (1 — np); 
the latter being a measure of the deviation of the Fermion 
occupation from the half-filled band situation and thus 
of hole doping. 

It is illustrative to compare this doping dependence of 
T*, derived from the dip in the DOS in the normal state, 
with the doping dependence of the mean field critical 
temperature T^ IF for amplitude fluctuation controlled 
superconductivity. That latter is characterized by the 
order parameters 

- = ^E<4 c *U p - ^T,(pt+pn (io) 

i i 

which refer to the off diagonal elements of the charge 
operators of the electron-pairs and hard-core Bosons re- 
spectively. Solving this mean field equation problem (for 
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FIG. 3: Evolution of the pseudo-gap near the chemical po- 
tential as a function of temperature for As = 1 and different 
doping concentrations n to t- 



FIG. 5: Doping dependence of T* , compared with the mean 
field critical temperature T C MF and the "phase fluctuation" 
temperatures T$. For comparison we also illustrate this dop- 
ing dependence by T^ B and T^ B when treating the Bosons 
as normal Bosons. 
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FIG. 4: Dependence of the minimum in the DOS on the con- 
centration of hard-core Bosons ub for different temperatures 
and Ab = 1. 



details for such an analysis the reader is referred to an 
earlier paper[Q) gives rise to a T C MF which exhibits a 
doping dependence which is quite similar to that of T* 
(see Fig. g), with an onset of amplitude fluctuation con- 
trolled superconductivity slightly below that tempera- 
ture where electron-pair fluctuations set in. T C MF has 
of course not the meaning of a transition temperature 
for the onset of superconductivity, which, as we shall see 
in the next section, is induced by phase rather than am- 
plitude fluctuations, except for the limit of low Boson 
concentrations where T C MF and T* smoothly join. 

It is an interesting question to ask how this doping de- 
pendence of T* changes when istead of hard-core Bosons 
one considers normal Bosons. We plot in Fig. ||the tem- 
perature T^ B , signalling the opening of the pseudogap in 
that case. It shows a monotonocally decreasing behavior 
with decreasing n tot , similar to that found for hard core 
Bosons but does not saturate, as is the case for those lat- 
ter, when approaching the fully symmetric limit n to t — 2. 
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V SPECTRAL PROPERTIES OF THE 
COOPERON PROPAGATOR 

Let us now examine the features of the normal state 
which act as a precursor of the superconducting phase. 
As previously shown J4fj, the intrinsically localized bound 
electron-pairs (Bosons), gradually acquire itinerancy as 
the temperature is lowered below T*. We shall here ex- 
plore this behavior as a function of doping and focus on 
the effect of the hard-core nature of those Bosons, which 
has been neglected in such previous studies. The res- 
onant electron-pair states in the Fermionic subsystem, 
induced by the exchange with the bound electron-pairs 
of the Bosonic subsystem, are described by the spectral 
properties of the Cooperon propagator 

C(q,T) = —Y, « c q-k T ( T ) c ^( T ) ; c k' l(0)cq-k' T (0) >, 

k,k' 

'11; 

which are intimately linked to the spectral properties of 
the single-particle Bose Green's function if(q, uj m ) via 
the relation 

C(q,w m ) = -^n(q,w m ) + -i-n 2 (q,o; m .)i ; i'(q, uj m ). (12) 

The thermodynamic and transport properties of our sys- 
tem are given by the low lying excitation spectrum of 
those resonant electron-pairs. Their spectral proper- 
ties are determined by spectral functions of the hard- 
core Bosons, given by the second term in Eq. 
In Figs. 



(12) 



we illustrate those spectral functions for 
the long wavelength regime together with their evolution 
with temperature for three representative concentrations, 
which cover the entire doping regime from underdoped 
(n to t — 1.65) with a high concentration of Bosons to the 
optimal/overdoped (n to t — 0.97) with a low concentra- 
tion of Bosons. We find that in the underdoped regime 
the Cooperons are well defined propagating modes with 
a narrow width of the spectral function which, moreover, 
strongly decreases with decreasing temperature (see Figs. 
||, 0). In the optimally/overdoped regime, on the con- 
trary, the spectral functions show overdamped mode be- 
havior (see Fig. ||). Tracing the peak position of the 
Boson spectral function as a function of wave-vector q 
permits us to determine the mass m p of those Cooper- 
ons. As we approach the dense limit of Bosons, m p in- 
creases sensibly, when we compare this mass for different 
values of doping (ntot) at a fixed given temperature (see 
Table I) . For the low doping regime where the Bosons are 
well defined quasi-particles, their DOS shows an evolu- 
tion with temperature in which the low energy part gets 
more and more peaked as the temperature is lowered and 
the peak position approaches the value Eo + bU(0, 0), (see 
Fig. ||, [l(]), as it should according to the Hugenholtz 
Pines theorempOj. 

We next turn to the evaluation of the concentration 
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FIG. 6: Temperature evolution of the Boson spectral func- 
tion for the low wave vectors for the underdoped regime with 
ritnt = 1.65. 
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FIG. 7: Temperature evolution of the Boson spectral func- 
tion for the low wave vectors for the underdoped regime with 
ntot = 1.20. 



of Cooperons, which, together with their mass will en- 
able us to estimate the phase stiffness in that system and 
thus the onset temperature of superconductivity due to 
phase fluctuations, Tq. The density of Cooperons, acting 
as superfluid charge carriers, is contained in the combi- 
nation of n p /m p entering the expression for the penetra- 
tion depth. Alternatively, and in an approximative fash- 
ion, it describes the density of itinerant quasi-particles 




FIG. 8: Temperature evolution of the Boson spectral function 
for the low wave vectors for the optimally/overdoped regime 
with That = 0-97. 



FIG. 9: Temperature evolution of the Bosonic DOS for the 
loi 



in the normal state, derivable from the Drude weight in 
the optical conductivity |5l|] . For the BFM scenario in- 
vestigated here, such a Drude component arises from an 
Aslamazov-Larkin term in the conductivity of the itin- 
erant electrons |40|| , involving the Cooperons and is con- 
tained in the second term of the Cooperon propagator, 
Eq. (12). In order to estimate the density n p of those 
Cooperons which give rise to such a Drude component, 
we have to attribute it to just that contribution of the 
Cooperon propagator, i.e., 

1 ^ 1 



N/3 



E 



(13) 



where the uncorrelated part 4j ^ q w IT(q, ui m ) has been 
subtracted out of the thermal average of the doubly oc- 
cupied sites, given by < c^c^cic^ >. In Table, II we 
present the mass and and concentration of Cooperons 
for different temperatures and doping rates correspond- 
ing to the well underdoped (n tot ~ 1.65) and the less well 
underdoped (ntot — 1.20) regime, for which the Bosons 
have well defined propagating quasi-particle features. 
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FIG. 10: Temperature evolution of the Bosonic DOS for the 
underdoped regime with ntot = 1.20. The insert presents the 
low energy part of this DOS. 



In order to highlight the effect of the hard-core na- 
ture of the Bosons we have repeated this study for the 
Cooperons for the case of normal Bosons rather than 
hard-core Bosons and present the corresponding values 
for the Cooperon mass and concentration by the numer- 
ical values, given in Table I, in the parenthesis. 



The above assessment of the mass and concentration 
of the Cooperons, contributing to the Drude peak in the 
normal state and ultimately to the superfluid current, 
shows that for hard-core Bosons their concentration de- 
pends little on doping, while their mass depends on it 
sensibly. The latter rapidly increases, as, upon reducing 
the hole doping, we enter the regime of high Boson con- 
centration where correlation effects become increasingly 
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TABLE I: Variation with doping n to t of m p (in units of jj) 
and n p for a fixed temperature T = 0.00667 and the estimated 
resulting T^. We present in parenthesis the corresponding 
values when the Bosons are treated as normal instead of hard- 
core Bosons. 



ntot 


m p 


n p 10 3 




1.651 


2.27 (1.29) 


9.60 (14.9) 


0.0042 (0.0116) 


1.614 


2.14 (1.24) 


9.49 (14.5) 


0.0044 (0.0117) 


1.495 


1.56 (1.06) 


9.04 (12.9) 


0.0058 (0.0112) 


1.376 


1.14 (0.88) 


8.44 (11.2) 


0.0074 (0.0127) 


1.204 


0.71 (0.62) 


7.30 (8.71) 


0.0102 (0.0140) 


1.156 


0.63 (0.55) 


6.91 (7.98) 


0.0111 (0.0145) 


1.105 


0.45 (0.49) 


5.88 (7.10) 


0.0124 (0.0145) 



important. As compared to the case of hard-core Bosons, 
for normal Bosons the variation with doping of the con- 
centration of Cooperons turns out to be much more im- 
portant, while for the mass it is less important. This 
leads, as we shall see below, to significant qualitative dif- 
ferences in the respective temperatures determining the 
onset of phase correlation driven superconductivity. 

On the basis of these findings, we now attempt to 
estimate the critical temperature for phase fluctuation 
controlled superconductivity. For that we put fcsT^ ~ 
Ti 2 (n p /m p )a, where n p and m p are our estimates for the 
density n s and mass m s of the superfluid charge carriers. 
a denotes a length scale which is of the order of the co- 
herence length or, alternatively, the inter-plane distance, 
depending on the degree of anisotropy of the system. 

Assuming a in the expression for to be given by the 
lattice constant, corresponding to a layer compound sys- 
tem such as the HTSC, we trace this critical temperature 
as a function of doping (ntot and 1 — Up) in Fig. 0. We 
notice a crossing of the energy scales related to the phase 
stiffness and to the electron pairing, as we approach the 
high doping limit, where the density of Bosons is small 
and Cooperons are no longer well defined quasi-particles. 
This phase diagram, Fig. ||, corresponds to that pro- 
posed on the basis of the phase fluctuation scenario Jl6[, 
but with the difference that there the doping dependent 
quantity of the phase stiffness was supposed to be re- 
lated to the density of superfluid carriers, while, accord- 
ing to our present findings, based on the BFM scenario, 
it should be primarily related to the mass of the super- 
fluid carriers - estimated as the mass of the Cooperons 
in our case. Upon approaching the optimal/overdoped 
regime, where crosses T* , the Cooperons loose their 
good quasi-particle features, and the onset of supercon- 
ductivity is becoming controlled by amplitude fluctua- 
tions, like in a BCS system. The opposite trend with 
doping of and T* , observed up to this level of doping, 
ceases accordingly and T c is constraint to decreases, since 
being limited from above by the decreasing behavior of 



TABLE II: Variation of n p , m p and tif with temperature T 
for ntot = 1-64 (top) and ntot = 1-20 (bottom). 



T 


n p 


m p 


n_F 


0.00667 


9.59 


2.23 10" 3 


0.995 


0.01000 


7.86 


2.63 10" 3 


0.994 


0.02000 


6.46 


5.81 10" 3 


0.990 


0.00667 


7.30 


0.71 10 -3 


0.988 


0.01000 


6.42 


1.00 10~ 3 


0.984 


0.02000 


5.56 


2.91 10~ 3 


0.971 



T* which controls pair formation. 

In comparison to these features, derived by considering 
hard-core Bosons, we find a noticeably different doping 
dependence of the phase fluctuation temperature when 
these hard-core effects are absent and the Bosons are 
treated as normal Bosons. See the corresponding values 
for T^ B in Table I and its graphical representation in 
Fig. [|. T^ B decreases only moderately upon approach- 
ing the low doping regime and, unlike for the case of hard 
core Bosons, does not vanish as the Boson concentration 
approaches ns = 1 for n to t = 2. 

In order to illustrate the combined effect of local elec- 
tron pairing and phase correlations of those resonant pair 
states, we now examine the Cooperon propagator, Eq. 
( |l2| ) in the low frequency limit and interpret it in terms 
of its physically intuitive form, given in eq. (2) in sec- 
tion II. Let us for that purpose illustrate this Cooperon 
propagator in Fig. [n] as a function of temperature for 
two cases, representing a well underdoped and a less well 
underdoped situations, showing the tendency with in- 
creased hole doping. Fitting the spatial dependence of 
the Cooperon propagator to this phenomenological form, 
permits us to determine the various parameter which 
characterize it and which we enumerate in Table III. We 
notice that upon increasing the temperature, the coeffi- 
cients C a , weighting the short range phase- uncorrelated 
local electron pairing, tend to rapidly decrease as we ap- 
proach T* in the entire doping regime. In the under- 
doped regime, upon decreasing the temperature and ap- 
proaching T^, the coefficients C^, weighting the phase 
correlated electron-pairs, increase rapidly together with 
the coherence length £, which is typically an order of 
magnitude bigger than the short range electron-pair cor- 
relation length rg. In the optimal/overdoped regime, 
on the contrary, hardly changes with temperature, 
while the coherence length follows the similar tempera- 
ture dependence as in the underdoped regime. A sys- 
tematic change is observed in the relative weight of the 
long range to short range contribution of the Cooperon 
correlation function: C^,/C a , which, with increased hole 
doping, shows a decreasing behavior when evaluated for 
some characteristic doping dependent temperatures such 
as T*. In particular, we find C^,/C a = 0.71, 0.45, 0.31 
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FIG. 11: Temperature evolution of the low frequency limit 
of the Cooperon propagator in the underdoped regime with 
Thtot = 1-65 and ntot = 1-20 (insert) and its spatial Fourier 
transform for n to t = 1.65. 



TABLE III: Characteristic parameters of the Cooperon prop- 
agator, Eq. (2), as a function of temperature T and for three 
doping rates n to t ~ 1-64 (top) and n to t = 1-20 (middle) and 
ntot = 1-02 (bottom). £ and ro are in units of the lattice 
constant a. 



T 


C a 




ro 


sen 


0.00667 


1.31 


1.06 


0.91 


13.75 


0.01000 


1.13 


0.61 


0.72 


7.59 


0.02000 


0.88 


0.42 


0.53 


3.58 


0.00667 


1.23 


0.66 


0.83 


11.76 


0.01000 


1.06 


0.45 


0.67 


6.90 


0.02000 


0.86 


0.37 


0.52 


3.51 


0.00667 


1.08 


0.33 


0.71 


9.07 


0.01000 


0.98 


0.32 


0.62 


6.18 


0.02000 


0.83 


0.32 


0.51 


3.45 



for ntot — 1-65, 1.20, 1.02 and corresponding values of 
T* ~ 0.016, 0.0142, 0.005 and T$ = 0.004, 0.010, 0.014. 
It is this relative weight increase of amplitude versus 
phase contributions, as we go from the underdoped to the 
optimal/overdoped regime, which indicates the change- 
over from phase correlation driven superconductivity to- 
ward amplitude correlation driven superconductivity. A 
remarkable result is that neither the short nor the long 
range scale depend sensitively on doping. We have some 
experimental indications |32"| from studies in the HTSC 
that, at least as far as £ is concerned, its doping depen- 
dence is very weak. 



VI EVIDENCE FOR PAIRING CORRELATIONS 
IN THE SPECIFIC HEAT 

The onset of the pseudogap, as seen in numerous exper- 
imental studies such as ARPES and single particle tun- 
neling, indicate a loss of low energy single particle spec- 
tral weight. This loss of single particle spectral weight 
ought to be accompanied by a compensating increase of 
spectral weight coming from collective excitations, which, 
for the present precursor scenario, should predominantly 
come from pair fluctuations. Without having to go to 
elaborate spectroscopic techniques, indications for such 
Many Body effects are already seen in basic thermody- 
namic quantities such as the specific heat CV(T) and 
entropy S(T), where a hump in Cy/T and a change in 
slope in S(T) is observed at temperatures around T*fl53|. 
A recent theoretical approa ch p4| on the basis of a classi- 
cal pair fluctuation scenario [55[ attributed this hump fea- 
ture very clearly to the contributions coming from pairing 
correlations, sitting on top of the single particle contri- 
butions. 

We shall in this section present a similar investigation 
on the basis of the two-component precursor scenario 
adopted here. We evaluate for that purpose the inner 
energy, given by 

U(T) = E[ ln (T) + Eg n (T) + Egf(T) 

k 

Eg n (T) = A B n B (T) 

i 

= -Jfa E n(q,w ro )if(q,a; m ) (14) 

and subsequently determine CV(T) = -^j,U(T) and 

S(T) = CdT'C v {T')/T'. The Fermion distribution 
function n^(T), the number of Bosons ub(T) and the 
expectation values of the interaction energy have to be 
calculated with respect to the full Hamiltonian. We il- 
lustrate in Fig. |l2] the temperature variation of Cy (T) 
and S(T) for different total doping rates n tot , correspond- 
ing to the underdoped situation and compare it with the 
non-interacting case (v = 0) for the case of very small 
concentrations of Bosons. The pseudogap is manifest in 
the dip-like feature of Cv(T) which occurs at T* , to- 
gether with a subsequent upturn of Cy(T) upon lower- 
ing the temperature which indicates the broad hump-like 
structure above T c . As we decrease ntot upon going from 
the underdoped to the optimal/overdoped regime, the 
dip moves to lower temperatures in agreement with a de- 
creasing T* and eventually disappears upon reaching a 
doping concentration where the number of Bosons tends 
to zero (given approximately by the non-interacting BFM 
with v = 0. At the same time the linear slope of Cy(T) 
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FIG. 12: Temperature evolution of the specific heat and en- 
tropy in the underdoped regime with ritot = 1.44, 1.20 and 
compared with the optimal/overdoped regime n to t = 1.02. 



at high temperatures increases in correspondence with an 
increase of the DOS at the Fermi energy and saturating 
at the value characterizing the low Boson concentration 
limit. Concerning the entropy, we are able to evaluate 
its increase, AS(T), starting form a fixed lowest temper- 
ature and going up to the highest temperature we have 
been considering. The rapid rise in AS(T) below T*, 
changing into a rather slow rise above T* signals the exis- 
tence of local order below T* which gradually disappears 
when going beyond T* to higher temperatures. 



VII SUMMARY 

A doping induced change-over from phase- to 
amplitude-fluctuation driven superconductivity is shown 
to result in a system with precursor pairing within a two- 
component scenario, involving charge carriers with differ- 
ent statistics: Fermions and Bosons, coupled together via 
a charge exchange term. The Fermions describe free elec- 
trons while the Bosons (more precisely hard-core Bosons) 
describe localized self-trapped electron-pairs having spin- 
| statistics which give rise to correlation effects in such 
a system. The opposite variation with doping of T* and 
T c is obtained, where the critical temperature is given by 
the phase stiffness of the system. With decreasing the 
concentration of the localized electron-pairs, the energy 
associated with this phase stiffness crosses the pairing 
energy ksT* in the itinerant electrons subsystem at a 
certain characteristic doping level. There, phase fluctua- 
tion controlled superconductivity changes over into am- 



plitude fluctuation controlled superconductivity, giving 
rise to a phase diagram, which, qualitatively, is reminis- 
cent of that proposed for the HTSC within the so-called 
phase fluctuation scenario. The doping dependent length 
scales for short range local electron-pair correlations and 
long range phase correlations are discussed on the ba- 
sis of the spectral properties of the Cooperon propaga- 
tor, describing the exchange induced pairing in the elec- 
tron subsystem. In the precursor pairing scenario studied 
here, it turns out that it is the degree of itinerancy of the 
Cooperons rather than their concentration which controls 
the doping dependence of the phase stiffness. This degree 
of itinerancy varies from well defined itinerant electron- 
pair states in the limit of high concentration of localized 
bound electron-pairs (low hole doping) to overdamped 
excitations in the limit of low concentration of localized 
bound electron-pairs (high hole doping). 

The phase diagram, Fig. ||, represented as a function of 
doping, involves changes in the concentration of electrons 
away from half filling (hole doping) for given changes in 
total concentration of charge carriers, ntot- The variation 
of the hole doping (changes in (1— np) of the order of 2%) 
in this phase diagram is small compared to changes in 
n t ot, which is due to the ID situation we have been con- 
sidering here. This result is less surprising when we con- 
sider a 2D system with an anisotropic charge exchange 
coupling between the bound electron-pairs and the bare 
itinerant electrons. One would then obtain a correspond- 
ing anisotropic pseudogap, similar to what is observed in 
the HTSC. Hole doping would now affect roughly equally 
all the regions near the 2D Fermi surface, and thus at- 
tributing only a very small fraction of the doped holes to 
the actual regions in the Brillouin zone where the pseu- 
dogaps are formed, i.e., around the M points and along 
lines parallel to [0,0] — [0,7r] and equivalents. 

One of the outstanding problems to be solved within 
such a precursor scenario for superconductivity is to un- 
derstand how the transition to the superconducting state 
occurs. This involves a competition between Cooper pair- 
ing and real space pairing and necessitates a generaliza- 
tion of the present system of Green's functions by in- 
cluding the anomalous Green's functions in order to treat 
superconducting fluctuations. This problem will be the 
issue of future studies. 
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